NASA Contractor Report 4226 


A Variable-Gain Output Feedback 
Control Design Methodology 


Nesim Halyo, Daniel D. Moerder, 

John R. Broussard, and Deborah B. Taylor 
Information & Control Systems, Incorporated 
Hampton, Virginia 


Prepared for 

Langley Research Center 

under Contract NAS1-17493 


NASA 

National Aeronautics and 
Space Administration 

Office of Management 

Scientific and Technical 
Information Division 


1989 



TABLE OF CONTENTS 


page 

LIST OF FIGURES v 

LIST OF TABLES vii 

I. INTRODUCTION 1 

CONVENTIONAL GAIN SCHEDULING 2 

II. VARIABLE-GAIN OUTPUT FEEDBACK-FORMULATION 6 

III. NECESSARY CONDITIONS AND INCREMENTAL COST 15 

IV. EMBEDDING INTO MULTI-CONFIGURATION CONTROL AND ALGORITHM 

DEVELOPMENT 23 

MULTI-CONFIGURATION CONTROL 23 

EMBEDDING INTO MCC 24 

ALGORITHM DEVELOPMENT 28 

V. APPLICATION TO RECONFIGURABLE AIRCRAFT CONTROL 32 

VI. CONCLUSIONS AND RECOMMENDATIONS 38 

REFERENCES 40 


PRECEDING PAGE BLANK NOT FILMED 




iii 


I 

I 


LIST OF FIGURES 


► 


page 


FIGURE 1. ELEVATOR INNER-LOOP CONTROL SYSTEM DIAGRAM 42 

FIGURE 2. EFFECTS OF THE GAIN SCHEDULED ON CLOSED-LOOP TYPE 0 

DFCS MAPPED EIGENVALUES 43 

FIGURE 3. CONTROL GAIN VARIATION WITH ANGLE OF ATTACK 44 

FIGURE 4. VARIABLE-GAIN OUTPUT FEEDBACK WITH DYNAMIC 

COMPENSATION 45 

FIGURE 5. AFTI F16 AIRCRAFT CONTROL SURFACE CONFIGURATION 46 

FIGURE 6. SIMULATION OF NOMINAL CLOSED-LOOP SYSTEM 47 

FIGURE 7. SIMULATION OF CENTERED LEFT HORIZONTAL TAIL AND 

RECONFIGURED CONTROL LAW 64 


PRECEDING PAGE BLANK NOT FILMED 

j g j i/ iMWMM'mm 


v 


LIST OF TABLES 


t 


i 


page 


TABLE 1. A AND B MATRICES AT 0.8 MACH AT 5000 ft\. ALTITUDE 81 

TABLE 2. CONTROL EFFECTOR TRANSFER FUNCTIONS 82 

TABLE 3. CONTROL EFFECTOR RATE AND POSITION LIMITS 83 

TABLE 4. CONTROL GAIN MATRICES 84 


WCtoTOi PACK H lank NOT FILMED 




vii 



I. INTRODUCTION 


In designing control laws, the usual first step is to describe the plant at a given 
operating point and then to develop a control law with a satisfactory performance for 
that plant model. In this case, an important property of a satisfactory control law is that 
its performance should not deteriorate in a material way when ‘small’ variations in the 
operating condition occur; i.e., the control law should be robust, or insensitive to small 
variations in the, usually unmeasured, parameters describing the operating condition. 

However, the design of most control systems of practical interest further requires 
control laws which maintain high performance in the presence of large changes in the 
operating point parameters. For example, the motion of an aircraft at a fixed airspeed 
can be described by well-known mathematical models [l]. However, the parameters of this 
model can vary considerably over a range of airspeeds. When dealing with large variations 
in the plant model parameters, a single constant-gain control law, no matter how robust, 
cannot attain the level of performance of a control law with variable gains. Thus, variable- 
gain control laws provide a class of controllers which can maintain high performance over 
a wide range of operating conditions. In Section II, we will formulate an optimal control 
problem for the design of variable-gain output feedback control laws.* 

Conventional gain scheduling techniques have provided a method of designing variable- 
gain control systems which can accommodate significant variations in the plant operating 
point parameters while continuing to make use of the accumulated knowledge and expe- 
rience in the design of linear systems. The ability to use the well-established theory and 

*The variable-gain output feedback problem was formulated in 1983: Halyo, N., “Modern 
Control System Design Using Optimal Gain Scheduling,” ICS TM-83-011, Information & 
Control Systems, Incorporated, 28 Research Dr., Hampton, VA, 1983. 
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accumulated experience in the design of linear systems, while extending its use and ap- 
plicability to control nonlinear systems makes the concept of variable-gain control laws 
highly attractive. However, the straight-forward application of the gain scheduling con- 
cept to complex systems can lead to unacceptable control system designs with deficiencies 
ranging from unsatisfactory performance to instability, as is discussed in the following. 

Conventional Gain Scheduling 

In its broadest form, the concept of gain scheduling can be applied to control systems 
designed using a classical frequency approach as well as a modern control approach. Mod- 
ern and classical designs for practical systems using specific simple forms of gain scheduling 
are not uncommon; e.g., [2], [3]. Figure 1 shows a simple schedule where only one gain 
is programmed as a function of one parameter. Currently, the use of any type of gain 
scheduling in modern control system designs is less common, [4], [5], probably due to the 
fact that few such systems are designed for practical applications, where the need for gain 
scheduling becomes apparent. 

The conventional approach to the design of control systems using gain scheduling 
consists of the following steps. 

1) A small number of operating points covering the operating range are selected; the 
method of selection is largely arbitrary, but the designer will use his intuition and 
previous experience in his choice. 

2) For each operating condition selected, a satisfactory linear design with the same basic 
structure is obtained using any available design technique. 

3) Using interpolation techniques, the control system gains are expressed as a linear 
function (i.e., a schedule) of the parameters describing the operating conditions, so 
that the schedule evaluated at any of the selected operating conditions is as “close” 
as possible to the design made in Step 2 for that operating condition. 

While the approach outlined above may result in satisfactory designs of variable- 
gain control systems (usually after a considerable amount of trial-and-error) , some of the 
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problems encountered in this approach are immediately apparent from the steps outlined. 
Some of these problems are given below: 

a) The selection of the points covering the operating range is quite arbitrary. Different 
selections can result in significantly different gain-scheduled control laws with differing 
performance characteristics. In general, the selected points are represented in the 
gain-schedule, while other points are not represented at all. 

b) The amount of effort involved in designing a satisfactory system for each operating 
point can be considerable. Particularly, if many points have to be selected to cover 
the operating range, this effort may become impractical. 

c) No matter what interpolation technique is used in Step 3, the gain schedule will not 
be the same as the system designed meticulously in Step 2 for the selected operating 
points, with the possible exception of unusual situations. Thus, even at the operating 
points selected, the gain-scheduled control system may not have a satisfactory per- 
formance. In fact, as demonstrated by the example below, the gain-scheduled system 
may even be unstable at one or more of the operating points selected. 

d) There is no guarantee that the conventional gain scheduling approach outlined will 
result in a satisfactory or even stable (!) gain-scheduled control law over the operating 
range even when such a system exists. 

e) The three basic steps outlined lack a control theory framework which provides insight 
into the solution of the problem resulting in an integrated variable-gain control law. In 
other words, each step is largely independent of the others; e.g., the systems designed 
for the particular operating conditions selected may be excellent designs for those 
conditions, but may collectively result in a very poor fit or an unacceptable gain- 
schedule. 

Conventional gain scheduling is essentially a curve fitting or interpolation problem 
rather than a control law design problem. The shortcomings of the current approach listed 
above are based not only on analytical observations, but also on the experience of the 
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authors in the design of gain-scheduled control laws for complex and practical applications. 
An example of some actual problems encountered can be found in [5] which considers a 
gain-schedule with respect to the parameters of angle of attack, airspeed and normal 
acceleration for a fighter aircraft. Figure 2 shows the closed-loop mapped eigenvalues of 
two gain-schedules obtained using a conventional gain scheduling approach. The individual 
operating point designs were obtained using modern control techniques. The interpolation 
method used selects the fit of the control gains which maximizes the correlation coefficient 
of the scheduled gains and the individual designs. From Figure 2, it is seen that Schedule 1 
results in an unsatisfactory, in fact unstable, design at one of the operating points selected, 
for the eigenvalues attributed to A/? command. Schedule 2, a simpler schedule in functional 
form, is seen to remain stable within the complete operating range. It is interesting to 
note, however, that the Schedule 1 gains are closer to the original design gains for the 
operating point considered than the Schedule 2 gains, as can be seen in Figure 3. Yet the 
closer gains are unacceptable from a stability standpoint. 

This example illustrates some of the remarks mentioned earlier. However, possibly a 
more important aspect is apparent here. The “closeness” of the schedule to the originally 
designed gains is not an appropriate criterion in the design of a variable-gain control system. 
Two gains may be “close” in most metrics (e.g., euclidean distance), but have very different 
stability characteristics, while gains which are not close can have more desirable and similar 
stability characteristics. 

Thus, the use of an arbitrary interpolation technique, based on minimizing a metric 
defined on the control gains, is seen to be inappropriate. A design methodology which, at 
least directly considers the stability of the gain-scheduled systems or, if possible, guarantees 
the stability of the scheduled system, would clearly be superior to the current heuristic 
approach. Furthermore, while stability is necessary for a control system to be acceptable, 
a satisfactory system must usually possess further desirable characteristics which depend 
on the particular control tasks to be performed. Thus, a design methodology for variable- 
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gain control system which addresses such desirable characteristics, at every point in the 
operating range, while guaranteeing that the resulting closed-loop system is stable in that 
operating range would eliminate the major problems discussed above. Such a control 
system design methodology will be proposed in the following sections in the context of 
modem digital control theory. 
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n. VARIABLE-GAIN OUTPUT FEEDBACK _ FORMULATION 


In this section, we will formulate an optimal control problem to design variable-gain 
output feedback control laws. The objective is to develop a control theory framework for 
the design of variable-gain controllers within which the shortcomings of conventional gain 
scheduling described in the preceding section can be overcome. In turn, this approach will 
result in a method of extending the operating range of the control law while continuing to 
use established linear control design and analysis techniques; i.e., a method of designing a 
controller for a nonlinear system using linear theory. 

Consider a nonlinear system linearized at a given operating condition defined by the 
parameter vector, p. 


x(t,p) = A{p) x{t,p) + B(p ) u(t,p) + w(t,p ) , (1) 

y{t,p) = C(p) x (t,p) + v(t,p) , (2) 

where p is a q-vector of parameters describing the system operating point, x(t,p) is a n- 
vector describing the plant state perturbation at time t for the given operating condition 
p, y{tfP) ^ a m-vector describing the measurements or the feedback variables at time t] 
w(t,p) and v(t,p) are vector random processes of appropriate dimensions describing the 
plant and measurement noises, respectively, at the operating condition p, and A(p), B(p), 
C(p) are matrices of appropriate dimensions which define the significant plant dynamics 
at the operating condition p. 

As we are interested in digital controllers, we will adopt a sampled-data control ap- 
proach to the problem. Thus, the control commands, u(f,p), will be assumed to remain 
constant over the sampling interval with a duration of T seconds. 
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u(t,p) = u(k,p) , kT <t < (k + l)T (3) 

The continuous stochastic system described by (1), (2) and (3) can now be described 
by the discrete stochastic system [6]. 

x(fc + l,p) = 4>{p)x{k,p) +T(p)u(fc,p) + w(k,p) , (4) 

y(k, p) = C(p) x(k, p) + v(k, p) , (5) 

where z(fc,p) is the continuous state at the sampling instant kT\ i.e., 

x(k,p)=x(kT,p) , (6) 

w(k,p) is the discrete plant noise sequence which is obtained from the continuous plant 
noise w(t,p) [6], and <f>{p), T(p), C(p) are the discretized versions of the continuous plant 
dynamics in (l) and (2). Note that if the control is of the form given in (3), then the 
discrete system in (4) describes the continuous state at the sampling instants, kT, with no 
approximation, as indicated by (6). 

When the continuous plant noise, w[t,p), is white, it can be shown that [6] the discrete 
plant noise is a sequence of uncorrelated random vectors; i.e., discrete white noise. In the 
following, we will assume that the discrete plant and measurement noises are both white 
sequences with zero mean; i.e., 


E[w(k,p)} = 0 , p)] = 0 

(7) 

E[w[k,p) w T {j,p)\ = W{p) 6 k -j 

(8) 

E[v{k,p) v T {j,p)} = V (p) 6 k -j 

(9) 
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where 6k is the Kronecker delta and the superscript “ T ” denotes the matrix transpose. 
Furthermore, it will be assumed that the plant and measurement noises and the initial 
state are uncorrelated; i.e., 


E[w[k,p)v T (k,p)] = 0 


( 10 ) 


E[w{k,p) x r (0,p)] = 0 , E[v(k,p) x T (0,p)] = 0 (ll) 

The control law structure of interest is the class of variable-gain output feedback 
controllers which may be described by 


= ~K{p)y{k,p) 


( 12 ) 


where y(k,p ) given by (5) represents the variables which have been selected for feedback by 
the designer and K[p) is the variable gain for the operating point, p. This controller struc- 
ture is a simple extension of the ‘constant gain’ output feedback case treated previously 


[7], [8]. 


In many cases, performance requirements make it desirable for the controller to use 
dynamic compensation in the feedback loop. Therefore, it is important to note that the 
output feedback control structure in (12) includes control laws with dynamic compensation; 
e.g., consider the controller structure 


z(k + l,p) = <t> z {p) z{k,p ) + T*(p) y(fc,p) 


(13) 


«(fc,p) = --K'z(p) z{k,p) - K y {p) y(k,p) (14) 

which contains the dynamic compensator (13) driven by the feedback vector y(k,p). The 
control law feeds back the compensator state z(k,p ) as well as the feedback vector y(k,p). 
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This controller structure is depicted in the functional block diagram in Figure 4. In this 
case, the matrices <f> z (p), r*(p), K z (p ) and K y (p) are the control law gains to be determined 
in the optimization of the controller. 

Select an open loop (pre-design) dynamic compensator 


*{k + l,p) = <t> zo {p) z(k,p ) + r, 0 (p) y{k,p) + u z (k,p) (15) 

where <f> zo {p) and r* 0 (p) are arbitrarily selected by the designer to achieve a desired ob- 
jective through the cost function. For example, <f> zo {p ) and r* 0 (p) can simply be selected 
as the appropriately dimensioned zero matrices. The vector u z (k,p) is the compensator 
control. Augmenting the state equations in (4) by the dynamic compensator in (15), the 
state vector x(k,p) by z(k,p), the feedback vector y(k,p) by z(k,p) and the control vector 
by u z (k,p), it can be shown that 


(u{k,p)\ = _( K u (p) K z (p) \ ( y{k,p) \ fl6l 

VMM/ V r Mp) -r,(p) <t> zo {p) -M p)) V*(M/ 

is the output feedback control law for the augmented system which produces the dynamic 
compensator control structure given by (13) and (14). Thus, this control structure is 
included in the output feedback structure of (12) by augmenting the system. 

Since some systems cannot be stabilized by instantaneous output feedback, the case 
of variable-gain dynamic compensator feedback is an important class of control laws. In 
the following, we will assume that the necessary augmentation to accommodate dynamic 
compensation has already been included in the system described by (4), (5) and (12). 

At a given operating point, p, the control structure (12) corresponds to the discrete 
stochastic output feedback formulation of [8]. However, in the current formulation, the 
quantity to be optimized is not a single control gain matrix for a particular operating 
point, but rather it is the global control law defined by the collection of control gains over 
the total operating range, {if(p),pe£}, where Z is the collection of all operating points, 


9 



p, of interest for the design under consideration. Therefore, the relationship of the control 
gains at different operating points must be an integral part of the optimal control problem 
we are formulating. 

Consider a linear functional relationship between the gain and the operating point 
parameters, i.e., 
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K{p) = K a + J2Pi K i . > ( 17 ) 

*=i 

where the coefficient matrices K 0 , • • • , K q determine the variable-gain control law over 

the total operating range. 

It should be noted that, since it is not feasible to compute the gain, K(p), at every 
point of Z, it is necessary to select some type of functional relationship between the gain 
and the parameters, whether linear or nonlinear. While the relation in (17) is linear in the 
operating point parameters, p,-, these parameters can be selected as nonlinear functions of 
the physical variables of the plant. Thus, the linearity constraint in (17) is considerably 
less restrictive than is at first apparent, but provides a structure which can be handled 
more easily in both analytical development and actual implementation. Nevertheless, a 
controller structure which can handle nonlinear relationships between the control gain and 
the parameter vector, p, more directly will be given below. 

The significant difference between the conventional gain scheduling approach and the 
variable-gain output feedback approach formulated here is that the latter approach does 
not design local control laws, but rather determines the global control law {K(p),pe Z} by 
optimizing a control objective over many operating points. If some operating points have 
less stability margin than others, the optimization automatically places greater emphasis 
on those operating conditions. 

Nonlinear relationships between the control gain and the operating condition param- 
eters can be treated in the following manner. Let Hi(p) be a r x r matrix of real functions 
of p, for each * = 1,2, As in (17), if, is a r x m matrix for * = 0,1, •••,?. Now 
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consider the relation 


K{p) = K 0 + Yl H T{p)Ki , peZ . (18a) 

»=i 

The elements of the matrices iT,(p) are arbitrarily selected integrable functions of p 
which may contain nonlinearities in terms of p. We will refer to this relation as a separable 
form, since the operating point parameter, p, and the gain coefficient matrices, K{, are 
separately expressed. The dual form 

<7 

K(p) = K 0 + Y,Ki G i{p) , peZ (18b) 

t=i 

where G,(p) isamxm matrix of real functions of p will also be considered. 

It should be noted that the linear relationship in (17) is a special case of the separable 
form in (18). The following selection of Hi(p) in (18) results in (17), 


Hi{p)=pil , 1 <t<q . (19) 

Thus, the variable-gain output feedback control law structure will be defined by (12) 
and (18) , resulting in a feedback law of the form 


( 20 ) 


u (fc.p) K. + Y,Hf{v)K i 

t=l 

Following the optimal quadratic sampled data formulation, consider a cost or objective 


y(fc,p) , pe£ 


function of the form 


Jic{K(p),p)= lim / E{x{t,p) T Q c (p)x(t t p) 

*/— o° tf Jo 

+ u(t,p) T R c {p)u(t,p)}dt (21) 

where E denotes the expectation operator. We will refer to Ji c as the local continuous 
objective function as it expresses the objective at the operating point, p. It can be shown 
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that the local continuous objective function Jt c is equivalent to a local discrete objective 
function [6]. In this study, we will use the form 

J t {K{p),p) = lim^ £ E{x{k + 1 ,p) T Q{p) x{k + 1 ,p) 

+ u(A:,p) T i2(p)u(fc,p)} (22) 

as the local objective function for the discrete optimization problem. For simplicity, cross- 
terms between state and control are not included in the objective functions (21) and (22). 

Since the parameter vector, p, has q components, the operating range, Z, is a subject 
of R q , where R denotes the reals. The operating range may be defined by the constraints 

< Pi < b{ , 1 < * < q (23) 

where p,- is the t th component of the parameter vector p, o* and &,• are the minimum and 
maximum values of that i th component, respectively. For notational convenience, we will 
express (23) by 


a < P < b (24) 

where in (20) “<” is a partial ordering on R q . 

In this formulation, the quantity to be optimized is not a single control gain for a 
particular operating point, but rather it is the global control law defined by the collection 
of control gains over the total operating range {K(p),pe Z}. Similarly, the objective or 
cost function for the design problem is not a local objective corresponding to the system 
performance at a single operating point, but a global objective with performance specifica- 
tions over the complete operating range, Z. Thus, the cost function for the design problem 
can be selected as 
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J=f *••/ f{p)Ji{K{p),p)dp q ---dp x - f f{p)Jt{K{p),p) dp (25) 

Ja i J a q Ja 

where f(p) is an non-negative scalar function of p selected by the designer to allow greater 
weighting of certain regions of the operating range over other regions. The shorthand 
notation on the right-hand-side of (25) uses a single integral sign from the vector a to b to 
denote q scalar integrals. 

In general, the operating range, £, need not be of the form of (23) and (24). Thus, in 
general, the desired global objective function for the variable-gain output feedback problem 
is 


J{K)= f f{p)Ji{K{p),p)dp (26) 

JJL 

where the local cost functions Ji(K(p),p) are weighted over the operating range, and 

(27) 

\ XVg / 

The notation on the left-hand-side of (26) simply recognizes that for the feedback structure 
selected, the global cost is determined by the choice of K . 

A discrete version of the global cost J can be expressed as 

M 

AK) = 'Ef> J ‘(K(P*hP i ) . ( 2 «) 

J = 1 

where Jt is defined by (22), fj is the discrete weighting for the j th operating point p 3 and 
{p 3 , j = 1,2, correspond to the plant operating points of interest for the design. 
This version of the cost function may also be considered an Euler approximation of the 
continuous cost function in (27). 

The variable-gain output feedback control law design can now be posed as a stochastic 
optimal control problem. The optimization consists of finding a variable-gain control law, 
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K(p), over the design operating range Z, which minimizes the cost function J or the 
discrete version J, subject to the constraints of (4) - (12) and (18). 

The optimal control approach taken here puts the variable-gain design problem into a 
firm theoretical setting. The desired characteristics of the controller are specified from the 
outset. No interpolation or curve fitting of local control laws arises in this approach. Since 
an instability at an arbitrary operating point would result in an infinite cost locally and 
globally for most practical systems, the optimal variable-gain control law will stabilize all 
the operating points considered in the design. Finally, the problem is formulated in a mod- 
ern control setting for direct-digital-design and is compatible with previous developments 
on stochastic output feedback [8], [9], [10], [11], [12], [13], [14]. 
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m. NECESSARY CONDITIONS AND INCREMENTAL COST 

In this section, we will obtain the necessary conditions for optimality of the variable- 
gain output feedback control problem posed in the previous section. Rather than using the 
Lagrangian approach and differentiating the augmented cost function to obtain the nec- 
essary conditions, we will follow the approach initiated in [8] of obtaining the incremental 
cost. From the latter, as in the stochastic output feedback and decentralized control cases, 
the necessary conditions will be apparent. 

At a given operating condition specified by the parameter vector, p, the system has 
the form of the standard stochastic output feedback problem defined by (4) - (12). 

Following [8], define the symmetric non-negative definite matrix P evaluated at the 
gain K[p ) as the solution of the discrete Lyapunov equation 

P(K(p)) =t T (K(p), p) P(K{p )) 4>(K(p),p) 

+C T (p)K T {p)R(p)K{p)C{p)+Q(p) , peZ (29) 

where 

4(K(p),p) = t(p) - r(p) K(p) C(p) , peZ (30) 

is the closed-loop transition matrix at the operating point, p. 

From Lemma 1 in [8], it is known that for all gains K(p) which stabilize the closed-loop 
plant, i.e., 

p(4,(K(p),p)) < 1 (31) 
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the Lyapunov equation (29) has a non-negative definite solution, P[K{p)).* In (31), p(<f>) 
denotes the spectral radius of the matrix, <f>\ i.e., the value of the largest magnitude among 
the eigenvalues of <f>. 

Furthermore, it can be shown that the local cost Ji(K(p),p) at the operating condi- 
tion, p, given in (22) can be expressed in terms of the non-negative definite matrix P(K(p)) 
when the closed-loop system is stabilized by the output feedback gain matrix, K{p). 

Lemma i. 

If p(<f>(K(p),p)) < 1, then the local cost Jt(K(p),p) defined in (22) is finite and is 
given by 


+ {k t W [r T M P(K(p)) r(p) + *(?)] k( p ) »( p )} (32) 

where P(K(p)) is the solution of (24). 

This results if a direct consequence of Lemma 2 in [8] and will not be proved here. It 
states that if the output feedback gain matrix K (p) stabilizes the open loop plant at p, then 
the limit in (22) converges to a finite cost Jt(K(p),p). When the plant (<£(p),r(p), C(p) 
is output stabilizable at p, then some gain, K(p ), stabilizes the plant and achieves a finite 
cost. If the plant is output stabilizable at all operating points of interest in Z, then some 
variable-gain {K(p),pe Z} will achieve a finite global cost for finite selections of Z. 


*Note that P depends both on the feedback gain K(p) as well as the operating condition 
p. For notational convenience, the dependence on p is not explicitly shown in the above, 
to avoid writing P(K{p),p). 
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We will define the stability sets 


S(p) = {K\p(t(K,p)) < X} , peZ (33) 

^{Wl.P^PW^W. VpeZ) (34) 

Thus, S (p) is the set of all output feedback gains matrices, K, which stabilize the plant 
at the operating point, p. Whereas S is the collection of variable-gain output feedback 
matrices which stabilize the plant at every operating point in Z. Clearly, the plant must 
be output stabilizable for all operating points of interest for 5 to be nonnull. 

To obtain the global cost functions given by (26) and (28), it is sufficient to integrate 
or sum the local cost in (32) over the appropriate points. Thus, 

W = \ J^f(p)ti{P(K(p))W(p)} dp 

+ IJ R m tr {* T (j«) P(K(p)) -K-(p) V(p ) } ip (35) 

JW = f £ U ftr {pW» Wtf)} 

3 = 1 L 

+ tr P(K(p>)) K{jP) V(p>')}] (36) 

where 

P(K(p)) = r T (p) P(if(p)) r(p) + R(p) , peZ . (37) 

Note that, following the notation used for P(K(p)), we neglect to explicitly show the 
dependence of P on p, but show the dependence on K(p), for notational convenience. 

In (35), the integration is over a subset of R q , namely over Z, and may be interpreted 
as q scalar integrals. In (36), the points {p*,j = 1,2, represent the operating 
points of special interest in Z. 
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At this point, the optimization of the variable-gain output feedback control is one 
of minimizing the global cost functions in (35) and (36) over the gains K which stabilize 
the closed-loop system. Since matrix Lyapunov equation solvers are readily available, the 
integrand in (35) and (36) can be computed at the desired values to perform the integration 
or summation. However, as the order of the systems considered increases, the numerical 
aspects of the optimization can be cumbersome. 

The incremental cost refers to the change in the cost due to a change in the control 
gains. For the local cost function Ji(K{p),p), the incremental cost at a given operating 
condition will be denoted by AJ(K(p), AK(p),p) and will be defined as 

AJ t (K(p),AK(p),p) = J,(K(p) + AK(p),p) - J t (K(p),p) , ptR , (38) 

where AK{p) is the change in the gain for the operating condition p. 

Let K{p) and K(p) + AK(p) belong to S (p). Then, the incremental cost can be shown 
to be 

AJ t (K(p),AK(p),p) = Itr (2Atf r (p) [p(/f(p) + AJf(p)) K(p)S(K(p)) 

- r T (p) P(K(p) + AK(p)) 4>(p) S(K(p)) C r (p)l 
+ A K T (p) P{K(p) + AK(p)) A K(p) S(K(p )) } , pe Z (39) 

where S(if(p)) is the solution to the Lyapunov equation. 

S(tf(p)) = <t>(K(p),p) S(K(p))<t> T (K(p),p) 

+ T (p) K(p) V(p) K T (p) r T (p) + W(p) , (40) 

5(Jf(p)) = C(p) S(K(p)) C T (p) + V(p) . (41) 
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As for P and P, S and S depend on K (p) and p, although the dependence to the latter is 
explicitly shown. It is interesting to note that S(K(p),p) is the steady-state covariance of 
the state x(k,p) when the control law K(p) closes the loop; i.e., 

S{K (p) , p) = lim E {x(fc, p) x T ( k, p) } . (42) 

To obtain the global incremental cost it is necessary to integrate (39) over the operat- 
ing points of interest. However, before taking that step, note that the variable-gain output 
feedback control structure which we have selected is of the form of (18). Now, let us form 
the matrix 


Then, 


H(P) = 


( B, (p) ) 


\H q {p)J 


peZ 


(43) 


K(p) = H r (p) K , 

peZ 

(44) 

AK(p) = H T (p) A K 

, peZ 

(45) 


Now, substituting (44) and (45) into (39), and combining the result with the global 
cost function in (35), we obtain 


AJ(K,AK) = J(K + AK) - J{K) 


(46) 


= itr |2 AJC t J f{p) H(p) P(K(p) + A K(p)) K(p)S(K(p)) 
- T T {p) P(K(p) + A K(p)) 4{p) S(K(p)) C T (p)l dp 
+ A K t J Up) H(p) P(K(p) + AK(p)) A K(p) S(K(p)) dpj 


(47) 
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Similarly, the incremental cost for the discrete cost function in (28) can be found to 
result in 


AJ(K,AK) = J(X + AX’) - J(X) 


(48) 


+ A K(f>)) K(pi)S(K(p<)) 

- r 7 V) P(K(p>) + A K(p>)) Mff) S(K(p>)) C T (p>) 

+AK T f^fj H(p>) P(K(p>) + Aif(p>)) AK(p 3 ) S(K(pi)) j (49) 

It should be noted that the expressions obtained in (48) and (49) are not approxi- 
mations such as first or second order variations of the cost, but rather represent the total 
change in the cost. Also note that the separable form selected for the gain in (18) is 
partially responsible for the form of the global incremental cost. 

From these expressions, the necessary conditions for optimality are easily obtained. 
Due to the form of the incremental cost, the gradient of the global cost due to a small 
change in the control gains, AX, can be obtained by observation. 


b r< 


M 


2AK T £ff(p>) 

y=i 


!£(<) = j z H(P) [?(*(P)) K(p) S(K(p)) 

- r T (p) P(K{ ,)) 4(p) S(K(P)) c r (p)| /(p) dp , (50) 

where K(p ) and A K(p) are given by (44) and (45). 

Setting the gradient to zero results in the necessary conditions. It may be of interest 
to partition the necessary conditions in the form: 
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[ H i (p)P(K(p))K(p)S(K(p))f(p)dp 

Jz 

= f Bi{p)T T {p)P{K{p))<f>{p)S{K{p))C T {p)f{p)dp 
Jz 

• *0, 1, — ,* (51) 


where 


H 0 [p) = / • (52) 

The necessary conditions for the linear case can be obtained simply by substituting 
(19) into (51). The gradient and necessary conditions for the discrete global cost function 
can be obtained following the same procedure explained above, resulting in 

§£«) = £ *&) \P{Ktf)) K(pi) $ 

- r T (p>) P(*y )) W) W)) cV)] /, (53) 


E J.V) £(*&>*)) W) SW)I /, 

y-t 

= E r T o») F(/f (p>)) 4 Kp') s(aV)) c V) ft , 

y-i 

* = 0,1, -.., 7 . (54) 

The solution of the necessary conditions would provide a potential solution to the 
optimal control problem, since any critical point of the global cost function will satisfy 
these conditions. Using the expressions developed for the gradient, it is possible to use 
standard gradient-based minimization techniques to obtain the optimal solution. 
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However, in the following section, we will show that the variable-gain output feedback 
problem can be embedded into the Multi-Configuration Control (MCC) problem which 
ICS has solved previously, and use the MCC algorithm to obtain the optimal gain. 
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IV. EMBEDDING INTO MULTI-CONFIGURATION CONTROL 
AND ALGORITHM DEVELOPMENT 


In this section, we will obtain a solution for the optimal control problem posed using 
the discrete cost function (28) by embedding it into an already solved problem, namely, the 
Multi-Configuration Control (MCC) problem [9], While standard minimization techniques 
can be applied using the expressions for the gradient obtained in the last section, these 
would not make use of any special knowledge about the form of the cost or incremental cost 
function. Whereas the approach used here makes use of some of the known characteristics 
of the cost function. 

From the familiar form of the expressions developed in the last section, it may be 
conjectured that the solution of the variable-gain output feedback problem may be simpli- 
fied. In fact, the similarities among the standard stochastic output feedback problem [8], 
the MCC problem [9], the decentralized control problem [9] and the variable-gain output 
feedback problem investigated here are largely due to the fact that the incremental cost 
function can be expressed in a similar form for each of these problems. 

Multi-Configuration Control 

First we will describe the Multi-Configuration Control (MCC) problem. The moti- 
vation for the MCC problem is the development of a modem control technique for the 
design of highly robust output feedback control systems; e.g., a single control law that can 
control a plant which has many configurations or many operating points. In comparison 
to the variable-gain controller, the MCC technique produces a constant-gain control law 
whose performance does not deteriorate as much as others when the plant operating points 
change without notice. Clearly, a variable-gain control law can provide better performance 
than a constant-gain control law when the operating point parameters are measured or 
estimated with sufficient accuracy. 
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Consider the plant described by (4) and (5). Suppose that we want to design a 
constant-gain output feedback control law which can operate not only at one nominal 
operating condition, but also at several other operating conditions. The form of the control 
law is given by 


u{k,p) = -Ky{k,p) , pc{pf t j = (55) 

where it is seen that the control gains do not depend on the operating condition. Therefore, 
the controller must perform satisfactorily, but at least be stable, at all the operating 
conditions considered. 

The optimal control problem posed by the global cost function given in (28) , the plant 
(4)-(ll) and the control structure (55), is referred to as the Multi-Configuration Control 
(MCC) or the Multiple Model Control problem [9]. The MCC problem has been treated 
and an algorithm to obtain the optimal MCC design has been presented in [9] and will not 
be repeated here. 

As can be seen from the formulations of the MCC and the variable-gain output feed- 
back control problems, the essential difference is in their control structures (55) and (12). 
Thus, to embed the latter into the former problem is a question of accommodating the 
control structures of these problems. 

Embedding into MCC 

The selection of the separable forms for the variable-gain control structure given in 
(18) provide the needed step for embedding. First, we will consider the separable form in 
(18a). For this case, recall that the gain can be written as 

K(p) = ff T (p) K , peR 
where H(p) is defined in (43) and K in (27). 
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Now, define the augmented control vector, u(fc,p), which is constrained to feed back 
only the variables y(k,p) through constant gains; i.e., 


*{k,p) = -K y{k,p) = - 


fK 0 \ 

K x 

\K a J 


y{k,p) 


(56) 


Note that the augmented control vector u(k,p) has (q + l)r components. 

In the same vein, define the augmented control effectiveness matrix, T(p), as 


l» = r( P ) h t ( p ) = (r(p) r(p)H?(p) ■ ■ • 


(57) 


With these definitions, note that 


f(p) U(k,p) = -r(p) H T {p) K y{k,p) 


(58) 


x(k + l,p) = <f>{p) x{k,p ) + T(p) u{k,p) + w(k,p) (59) 

Thus, the plant state remains the same, while the control vector is augmented but 
still restricted to the feedback vector, y(k,p), given by (5). 

Finally, define the augmented control cost matrix, R(p), by 


m = mp) r(p) h t (p) = 


( m 

n,(p)R(p) 


R(p)HT(p) 

Hi(p)R(p)hT(p) 


\H,(p)R(p) 


R(p)Hf(p) \ 

Hq(p)R(p)Hj(p) y 
( 60 ) 


It follows that the local, hence global, cost functions in (22) and (28), can be expressed in 
terms of the augmented control vector in a quadratic form. 
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( 61 ) 


1 N 

MK(p),p) = foa* 2 , N + ^ 53 E i x ( k + l 'P) T G(p) x ( k + 1 »p) 

' ' Jk=o 

+ u(fc,p) r J R(p)u(fc,p)} 

j(/o = £/y WpO.pO ( 62 ) 

j = i 

The augmented plant model (59), the feedback vector (5), the constant-gain control 
law of (56) and the global cost in (62) define a Multi- Configuration Control design problem, 
which can be solved using the MCC algorithm. The optimal control gain matrix, say K* , 
obtained for this problem also provides the optimal variable-gain control law, say K*(p), 
for the problem posed in Section II through the relation 

=H T {p)K- (63) 

The dual of the previous development is obtained when the separable form in (18b) 
is used as the control structure. In this case, rather than augmenting the control vector, 
the feedback vector is augmented. Thus, let 


•••if,) , 


(64) 


G(p) = 


( 1 \ 
<?.(?) 

W(P)' 


K(p) = K G(p) 

Define the augmented feedback vector, y(k,p), as 


(65) 


( 66 ) 


J(lc,p) = G(p) y(k,p) = C(p)x(fc,p) + v(fc,p) 


(67) 
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<?(p) = G(p) C(p) = 


( C(P) 

capjcw 


VC,(p)C(p)/ 


0(k,p) — G(p) v(k,p) 


( 68 ) 


(69a) 


V(p) = P) $ T (k,p)} 


(69b) 


K(p) = G(p) V[p) G T (p) = 


/ ^(P) 

Gi(p)V(p) 


^(p)Gf(p) 


VG,(p)V(p) 


^(PJG^P) \ 

GMV(p)Gf(p)J 


(70) 


u(b.p) = -£y{k,p) ■ (71) 

The plant model given by (4), the feedback vector given by (67)- (70), the control 
constraint given by (71) and the cost function in (28) pose an optimal Multi- Configuration 
Control problem whose solution, say K* , provides the optimal variable-gain output feed- 
back control law in the form of 

k’{p) = Z‘g(p)=k; + Y, k; g, (p) • (72) 

»=1 

The optimal variable-gain output feedback problem has thus been embedded in the 
MCC problem whose solution can be obtained by the algorithm described in the following 
when the original problem is augmented as described above. It should be noted that 
this development has been obtained at the expense of augmenting the control or feedback 
vectors according to which form is used, which increases the dimension of the corresponding 
vector. This increase in dimension is proportional to the number of elements used in the 
separable form (18) of the control law. Thus, the complexity of the selected variable-gain 
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control law proportionally determines the dimension, hence the numerical complexity, of 
the problem to be solved. 

Algorithm Development 

As described in the preceding, the variable-gain output feedback problem must first be 
transformed into a Multi-Configuration Control (MCC) problem. Then, the MCC output 
feedback algorithm can be used to obtain a critical point of the global cost function. The 
following algorithm can be used to obtain the gain K* defined by (64). 

Variable-Gain Output Feedback Algorithm 

1. Embed the variable-gain problem into the MCC form by augmentation, using the 
augmented constant gain K. 

2. Select an initial stable gain R 0 , oc 0 = 1 ,z> l,t = 0. 

3. Solve the Lyapunov equations j = 1, 2, • • • , M. 

4> ! {R i ) = <t>(£<G(p)p’) , Tj = IV) , Cj = c(pi) 


Q, = Q(p’) , Rj = B(p’) , w, = v*V') . v# = «V) 


Pitfi) = M&) T P,i&) M&) + cjRfRj K Cj + q, 


Sj(R i) = tjfRjSjtfi) 4,(Ri) T + r, RiVj Rf rf + tv,- 


If Pj{£i) or Sj(Ki) is not non-negative definite, go to 6. 

4. Solve for the direction d(Ki) 


M 
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a QJ 

where Pj(Ki), Sj (£,) and -=-(£,) are given in (37), (41) and (53), respectively and 

dK 

fj is the discrete weighting of the } th operating point in (28); now compute the new 
gain 


Ri + 1 = Ri + a. d(j£i) 

5. Compute the cost J(/C»). If * = 0, set i = 1, go to 3. If 


then go to 7.* 

6. Reduce step 


1) <o 


«i = ajz , (£,-) = (£,_,) , <J(r<) = <!(«(_ ,) 


(£,-i) = (£,) +“i <*(£.) , a i+1 = a, , i = i + l , go to 3. 

7. Check convergence criteria; if not converged, go to 3. 

The algorithm described above was coded for the case of linear relationship between 
gain and operating conditions; namely, <?,(/>) = p{, i = 1,2 ,*•*,?• When the number 
of controls and the number of feedback variables in the problem are large, solving for the 
iterative direction, d(£»), in step 4 is quite difficult and computationally costly. In variable- 
gain output feedback problems, these dimensions are often large due to the augmentation 
in transforming the problem into the MCC form. 

*This condition may be replaced by the alternative condition 

J( Jfi) - J ( < a - l(2 ~ ni ~ l) E", fj tr , 

which would produce a greater improvement per iteration in the absence of numerical 

errors. 
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To increase the speed and accuracy of the algorithm, a new method of solving the 
direction d(K) was developed. First, note that the equation can be rewritten in terms of 
Kronecker products, ®: 

Y u (R) ® A(^) 

>=i 

where v{»} is the column vector form of the gain matrix in brackets. Thus, obtaining d(£) 
requires the inversion of a large matrix in brackets in (73). To reduce the computational 
load involved, first we will make the following approximation by minimizing 

M 

|| ®P, -S®F|| 2 , (74) 

/» 1 

where the argument K has been dropped for convenience. When S and P are known, the 
required inverse is easily obtained by noting that 

(s®p)~ l =5" 1 ®P" 1 (75) 

After considerable manipulation, it can be shown that S and P satisfy 

1 M 

5 E s> py ( 76 > 

i M 

P = -4— Y fj Pj tr{& S } (77) 

Numerical experience indicates that iterating on (76) and (77) results in robust and 
rapid convergence to a solution which is a sufficiently good approximation. However, even 
with this approach, inverting S and P can be difficult when S is obtained by augmentation 
from the variable-gain problem. S is then dimensioned ( q + l)m x (q + l)m. To alleviate 
this problem, we minimize 


v {<*(*■)} = -v 1 1^(£) } (73) 
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I|5 - 5, ® Sjf 


(78) 


where Si is dimensioned (g -f 1) X (g + 1) and S 2 is dimensioned m x m. 

We will use the following notation. The i,j element of any matrix A will be denoted 
by [A],-y. The (g + l)m x (g + l)m matrix 5 will be partitioned as shown below: 


S = 


( 5(1,1) 5(1,2) 

5 ( 2 , 1 ) 


S(x,«+1)\ 


J 


(79a) 


V 5(g + 1, 1) 

For given i and j, S(i,j) is a m x m matrix; so that {S(t,j), 1 < t, j < g+ 1} determines 
all the elements of S. In the same vein, for given k and £, define T(k, l) as the (g+1) x (g+1) 
matrix given by 


[T(A:,£)] t y = [5(i,j)]^ , 1 < *', J < g + 1 , l<k, t<m (79b) 

Now iterating alternately on (80) and (81) results in 5j and S 2 which provides an 
approximation to 5. 


= TT^-jT^ tr (52 5(*,i)} 

lP2ir 


i < *, y < 9 + 1 , 


(80) 


IPill 


(81) 


It should be noted that 5 _1 can be obtained by directly inverting 5 whenever desired 
or by using the approximation described above coupled with 


S~ l cz 5f 1 ® &2 1 

Withe these improvements, the impact of high dimensionality can be significantly 
alleviated, making even high-order variable-gain output feedback design problems solvable 
with reasonable effort. 


31 


V. APPLICATION TO RECONFIGURABLE AIRCRAFT CONTROL 


The design methodology and algorithm developed in the preceding sections was ap- 
plied to a high-performance aircraft reconfiguration problem arising from a control surface 
failure. The failure of a control surface, a control actuator or servo, or the failure of 
a sensor in the feedback loop provide excellent examples to illustrate the control design 
methodology developed, although numerous other applications to control system design 
are possible. Restructuring and reconfiguring aircraft flight control law has received some 
attention [15], [16], [17]. 

The high-performance aircraft considered in this application is the AFTI F16 shown 
in Figure 5. The objective is to design a reconfigurable digital control law in which a failure 
in the horizontal tail is accommodated by the automatic reconfiguration of the flight control 
law as soon as the failure is detected. While the detection of the failure is an integral part 
of the process, here we concentrate on the design methodology for the control law. 

As the AFTI F16 has static instability, the design of the control law and its re- 
configuration strategy gain more importance, since a failure can easily bring the aircraft 
back to instability with catastrophic consequences. The objective of designing a digital 
reconfigurable controller will be attempted using the variable-gain output feedback design 
methodology developed in this investigation. 

The AFTI F16 aircraft rigid-body dynamics and servoactuators were simulated on a 
digital computer. The reconfigurable control system design was formulated as an opti- 
mal variable-gain output feedback control problem. A reconfigurable design was obtained 
using the algorithm developed in the last section. Finally, the closed-loop reconfigurable 
aircraft/controller system was simulated to evaluate its behavior. 

The simulation of the aircraft and its servoactuators was done using a linearization 
about a Mach number of 0.8 and an altitude of 5000 ft. with the aircraft in nominal 
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straight and level flight. At this flight condition, the aircraft has a static instability in the 
short period mode at +0.3172. Furthermore, since the aircraft uses a single engine, the 
gyroscopic effects are not negligible and produce a roll/yaw coupling. As a result, even the 
case of no control failure contains significant coupling between the lateral and longitudinal 
modes. These effects are readily apparent in the coupled (A, B) matrix model shown in 
Table 1. These A and B matrices were obtained by linearizing the NASA Langley Research 
Center AFTI F16 nonlinear simulation. 

The aircraft control surfaces can be seen in Figure 5. The independently movable 
controls consist of the rudder, left horizontal tail, right horizontal tail, left and right 
flaperons, left and right vertical canards and the throttle position. The leading-edge flaps 
and the speed brake were not used in this simulation. Detailed integrated servoactuator 
models were used to simulate the dynamics of the control surfaces. A first order model 
with a 0.5 sec. time constant was used to model the throttle command to actual throttle 
position. The remaining model transfer functions are shown in Table 2. The control 
effector simulation uses both rate and maximum/minimum position limits which are shown 
in Table 3. 

The simulation uses a 3 rd -order Adams-Bashforth numerical integration algorithm to 
update the aircraft variables. The simulation update or sampling rate is 40 Hz. The 
complete rigid-body dynamics and servoactuator simulation model is of 33 rd order. 

To design a variable-gain output feedback control law for the aircraft, first requires 
building a design model. The design model determines the structure of the control law as 
well as the values of the control gains which are used. We start building the design model 
with the aircraft rigid-body state and control variables. 

x T = (ub tiifi qs 0 vb Pb r B <t> VO (82) 

where ub,vb,wb are the linear velocity components along body xb , Vb > Z B axes, qBtPB^B 
are the pitch, roll and yaw rates expressed in the body axes, and are the pitch, roll 
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and yaw angles. 

The control effector position vector will be denoted by 6 with the following order of 
variables. 


6 t = {6R Shtl Shtr Sfl Sfr Svel Sver 6th ) (83) 

where SR denotes the rudder position, Shtl and Shtr the left and right horizontal tails, 
Sfl and 6 fr the left and right flaperons, Svel and Sver the left and right vertical canards, 
and 6th the throttle position. For a variety of reasons (e.g., see [13], [10]), we select a 
control rate command structure for the controller by augmenting the design model by the 
following equations. 


6k+i — &k + ujfc > (84) 

where u* is the control vector for the design model and represents the change in the 
corresponding control effector position from one sample to the next. A sampling rate of 
10 Hz is used for the purpose of obtaining the control law design. 

Finally, to obtain a reconfigurable design with type-1 steady-state characteristics, 
integral feedback of the command variables is added. The command state vector, z is 
selected as 


z T = (uj s 6 <f> ip) , (85a) 

Z T = (“Be 4>c V»c) (85b) 

where the subscript “ c ” denotes the commanded value of the variable. The integral error, 
or more correctly accumulated error, is obtained by 

h+i =h + A t{z k - z ck ) (86) 
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The complete design model state, X, is the augmented state shown below 


X T = {x T S T I t ) , (87) 

with the control vector being denoted by u. The discretized equations for the aircraft 
dynamics, augmented by the (84) and (86) forms the design model state equations. 

f Xk+i\ ( <i> T(p) (A / x k \ [0\ ( 0 \ 

X k+1 = \ 6 k+l = 0 / 0 U* + J u* - 0 (88) 

\h + ij \AtH, 0 lj\hj \0j \Atz ck J 

r(p) = r 0 -r 1 p 1 -r 2 P2 , o< Pt <i , * = 1,2 ( 89 ) 

In the equations above, the parameter vector p has two components, with pi and p 2 
representing the fractional surface effectiveness loss for the left and right horizontal tail 
surfaces. In (89), T 0 is the standard matrix obtained by discretizing the continuous aircraft 
perturbation equations when all the controls are operating normally. Ti is a null matrix 
with the exception that its 2 nd column (i.e., the column corresponding to the left horizontal 
tail) has been replaced by the 2 nd column of T 0 . Similarly, T 2 is a null matrix with the 
exception that its 3 rd column (i.e., the column corresponding to the right horizontal tail) 
has been replaced by the 3 rd column of r o . Thus, p, = 0 corresponds to zero effectiveness 
loss or a normally operating control surface; whereas a value of p, = 1 corresponds to 100% 
loss of effectiveness such as a centered surface. 

The variable-gain control law is of the form 


K{p ) = K 0 + K\pi + K 2 P 2 

( 90 ) 

u k = —K(p) X k 

( 91 ) 
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Figures 6 and 7 show the simulations of the reconfigurable control law obtained with 
the design model described and using the conditions where 1) all controls working normally 
(pi = P 2 — 0), 2) the left horizontal tail is centered (pj = 1 ,p 2 = 0), and 3) the right 
horizontal tail is centered (p i = 0,p 2 = 1). Table 4 shows the gain matrices used in the 
simulations. 

The command in these simulations is a total pitch angle of 15°, while the commanded 
roll and yaw angles are zero, and the commanded speed is Mach .8. In Figure 6, the case 
where no failures occur is simulated. Figure 7 shows the simulation for the case where the 
left horizontal tail is centered at 1.8 sec. into the simulation while the commanded maneuver 
is in progress. The simulation assumes that the failure is immediately detected and isolated; 
so that the parameter vector, p, immediately reconfigures the flight control law. While a 
more realistic FDI simulation would be desirable, as long as the FDI correctly isolates the 
failure, the difference from the results shown here would be in the transient behavior while 
the steady-state results would be the same [18], at least for linear simulations. 

In Figure 6, with no failures occurring, the aircraft pitch angle rapidly increases from 
its trim value at straight and level flight, slightly overshoots the commanded 15° and then 
smoothly settles to its commanded value. 

The maneuver is largely achieved by the horizontal tail producing a positive pitching 
moment. The pitch-up results in more drag and a correspondingly small drop in speed 
which is then controlled by an increase in the thrust. However, the coupling between the 
thrust and the rolling and yawing moments produces a slight response in the lateral modes. 
This coupling produces a maximum of 0.09° of roll and 0.03° of yaw with an oscillatory 
behavior which is quickly damped out. Further experimentation with the design would be 
needed to eliminate this transient effect. 

When during the initial stages of the commanded maneuver, at 1.8 sec., the left 
horizontal tail fails and is automatically centered, the variable-gain controller has to re- 
configure its strategy in order to accommodate this condition. Despite this failure, the 
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reconfigured control law easily meets its objective of 15° pitch with no steady-state roll or 
yaw. Although the response in pitch is different, it is just as fast as before. The roll and 
yaw angles go through a larger transient, but are firmly brought back to their commanded 
values after the reconfiguration in the control law. 

The failure of the left horizontal tail (but not the right) produces a large coupling 
between the lateral and longitudinal dynamics. The control law reconfiguration strategy is 
to use the operational right horizontal tail to produce and maintain the pitching moment 
necessary to achieve the commanded pitch, while nulling the cross-coupling moments with 
the flaperons and vertical canards. The thrust profile seems only slightly changed. 

After the failure, the large deflection in the right horizontal tail, now unbalanced, 
produces a significant amount of positive rolling moment as well as positive yawing mo- 
ment. The sudden onset of these moments, produce the positive roll and yaw transients 
mentioned earlier. The reconfigured control law is seen to bring and maintain the roll at 
zero mainly by using the flaperons to counterart the rolling moment, and by using the 
rudder to counteract the yawing moment. The vertical canards are used to shape and 
improve the transient behavior. 

Overall, the reconfigured control strategy is seen to handle the failure of the left 
horizontal tail with relative ease. Alternately, the design methodology of optimizing a 
variable-gain output feedback controller has produced a reconfigurable control design which 
accommodates the control surface failure considered with relative ease. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


The main contribution of the investigation described in this report has been the for- 
mulation, development and solution of the variable-gain output feedback problem in the 
form of an optimal stochastic control problem. This approach provides a control the- 
ory framework within which the operating range of a control law can be significantly 
extended. Furthermore, the approach avoids the major shortcomings of the conventional 
gain-scheduling techniques. 

The optimal variable-gain output feedback control problem is solved by embedding 
into the Multi- Configuration Control (MCC) problem, previously solved at ICS. An al- 
gorithm to compute the optimal variable-gain output feedback control gain matrices has 
been developed. The algorithm is a modified version of the MCC algorithm improved 
so as to handle the large dimensionality which arises particularly in variable-gain control 
problems. 

The design methodology developed was applied to a reconfigurable aircraft control 
problem. A variable-gain output feedback control problem was formulated to design a 
flight control law for an AFTI F16 aircraft which can automatically reconfigure its control 
strategy to accommodate failures in the horizontal tail control surface. Simulations of the 
closed-loop reconfigurable system show that the approach produces a control design which 
can accommodate such failures with relative ease. 

While the example considered is an important illustration of the power of this new 
design methodology, applications to a large variety of current problems is desirable. In 
particular, it is possible to extend the flight regime of most aircraft by appropriate appli- 
cation of the methodology. For example, superagility characteristics can be achieved by 
using angle-of-attack and airspeed as components of the operating condition parameter, 
p. Sensor failure accommodation can also be achieved using the methodology. Numerous 
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other applications to a variety of control problems remain for future investigation. Two 
important areas for study are: 1) the development of new algorithms that have better 
numerical convergence characteristics for high dimensional problems, and 2) the extension 
of the variable-gain output feedback approach to the feedforward control law design. 

Finally, it should be noted that digital control design methodologies for nonlinear sys- 
tems are rare at present. Until such nonlinear techniques and analysis tools able to handle 
the theoretical and, more importantly, practical requirements of complex control systems 
can be developed and demonstrated, the variable-gain output feedback approach presented 
here provides the control system engineer with a design methodology for aerospace control 
problems. 
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Figure 1. Elevator Inner-loop Control System Diagram 
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Figure 4. Variable-gain Output Feedback with Dynamic Compensation 
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Figure 6.a. Simulation of Nominal (No Failure) Closed-loop System 




880.00 



1 1 1 1 1 1 1 1 1 1 
o 

i i i i i | i n i i 
O 

i m ii 1 1 i i i i 

o 

i 1 1 i T'| i r m 
O 

1 1 1 up 

o 

o 

O 

o 

O 

o 

to 

d 

in 

o 

to 


h- 

CO 

to 

to 

00 

00 

00 

00 

00 


(oas/ij) loian 


o 

Fo 


-O 

-o 

-00 


= o O 
-o 0 ) 
-<6 






E-< 


-o 

-<D 
- <N 



50 


Figure 6.d. Simulation of Nominal (No Failure) Closed-loop System 
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Figure 6.e. Simulation 
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Figure 6.1. Simulation of Nominal (No Failure) Closed-loop System 
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Figure 6.m. Simulation of Nominal (No Failure) Closed-loop System 





Figure 6.n. Simulation of Nominal (No Failure) Closed-loop System 
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Figure 7.b. Simulation of Centered Left H< 
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Figure 7.d. Simulation of Centered Left Horizontal Tail and Reconfigured Control Law 
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7.m. Simulation of Centered 
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Figure 7.n. Simulation of Centered Left Horizontal Tail and Reconfigured Control Law 
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Figure 7.o. Simulation of 




Table 1. A and B Matrices at 0.8 Mach and 5000 ft. Altitude 
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A. Integrated Servoactuator (ISA) Transfer Function 

G(s) = (20.2) (144.8) (71.4) 2 
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Rudder 

Horizontal Tails 
Flaperons 
Vertical Canards 


B. Leading-Edge Flaps Transfer Function 
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s + 7.38 



C. Speed Brake Transfer Function 


<3 M = 


D. Percent Throttle Transfer Function 
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G(s) = 
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Table 3. Control Effector Rate and Position Limits 


Control 

Rate Limit 

Max. Positive 
Deflection 

Max. Negative 
Deflection 

Rudder 

120 deg/s 

30 deg 

30 deg 

Left Hoizontal Tail 

60 deg/s 

25 deg 

25 deg 

Right Hoizontal Tail 

60 deg/s 

25 deg 

25 deg 

Left Flaperon 

52 deg/s 

20 deg 

23 deg 

Right Flaperon 

52 deg/s 

20 deg 

23 deg 

Leading Edge Flaps 

30 deg/s 

25 deg 

2 deg 

Left Vertical Canard 

108 deg/s 

27 deg 

27 deg 

Right Vertical Canard 

108 deg/s 

27 deg 

27 deg 

Speed Break 

27.28 deg/s 

60 deg 

0 deg 

PLA 

N/A 

100% 

0% 




Table 4. a. Control Gain Matrices 
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